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19.  Abstract  (ODntinued) 

In  Reference  1,  it  was  postulated  that  LDP  was  also  responsible  for  phase  changes,  in¬ 
creased  chemical  reactivity,  and  anomalous  electrical  activity,  which  are  experimentally  ob¬ 
served  in  certain  inert  materials  under  relatively  low-level  shock-loading.  This  postulation  is 
supported  by  the  present  LDP  assessment  results  for  shocked  PPMI,  iron,  and  NaCl.  As  such, 
Fitzgerald’s  LDP  Vf  is  strongly  believed  to  be  the  “trigger  mechanism”  responsible  for  the  low' 
pressure  “catastrophic  shock”  and  consequent  “bond  scission”  event  suggested  by  Graham  [26], 
which  was  based  on  many  experimental  observations. 


The  importance  of  the  De  Broglie  momentum-wave  length  relation  with  respect  to  particle 
vclccitv  propagation  in  shock  IcadoJ  Suliu  maitii-Is  is  emphasized  via  the  results  reported  in 
this  document. 
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1.  INTRODUCTION 


Reference  1  documents  the  application  and  extension  of  Fitzgerald’s  crystal  lattice  disinte¬ 
gration  criteria  [2,  3,  4]  to  shock-loaded  solid  reactive  and  inert  materials  In  Reference  1,  six 
explosives  (comp-B3,  comp-B,  pressed  TNT,  tetryl^H-6,  and  PBX-9404)  and  one  inert  material 
(PMMA  or  plexiglas  II)  were  investigated.  It  was  tentatively  concluded  that  crystal  lattice  disin¬ 
tegration  was  responsible  for  the  experimentally  observed  response  which  occurred  at  certain 
shock  pressure  (Ps)  levels. 

To  further  investigate  this  initial  conclusion,  additional  available  experimental  data  for 
several  reactive  (explosive)  and  inert  materials  have  been  examined,  and  the  important  particle 
velocities  (Vj  and  Vf)  were  computed.  The  associated  shock  velocity  (Usf)  and  shock  pressure 
(Psf)  were  also  calculated  to  facilitate  comparison  with  experimentally  observed  reactive  or 
anomalous  events. 

The  explosives  analyzed  are:  pressed  TNT,  PBX-9407,  octol,  pressed  PETN,  and  crystal 
PETN.  Pressed  TNT  was  considered  again  because  some  important  experimental  information 
was  found  since  the  publication  of  Reference  1. 

Three  rather  diverse  inert  materials  are  considered.  They  are:  poly  pyromellitimide 
(PPMI),  iron,  and  sodium  chloride  (NaCl).  PPMI  is  a  complex  polymeric  compound  whose 
low-pressure  shock-loading  behavior  is  highly  non-linear.  Iron  and  NaCl  have  both  been 
extensively  investigated  via  shock-loading  experimentation.  The  results  have  been  controver¬ 
sial,  unexpected,  and  usually  inexplicable,  since  both  of  these  rather  common  substances  exhibit 
exceptionally  complex  behavior  under  shocked  conditions. 

The  present  suggestion  that  Lattice  Disintegration  Phenomena  (LDP)  or  atomic  bond  break¬ 
ing  Vf  effects  are  a  “root  cause”  for  low  pressure  shock  induced  reactions  may  also  be  controver¬ 
sial.  However,  the  number  of  favorable  comparisons  of  predicted  and  experimental  information, 
contained  in  thic  repcrt  :.r.d  in  Rcfercr.wt  1,  is  mote  ihuu  could  reasonably  be  expected  from 
mere  chance  or  fortuitous  circumstances.  Four  of  the  nine  different  explosives  have  been  ana¬ 
lyzed  for  two  different  densities,  which  makes  a  total  of  13  comparative  cases  for  the  energetic 
materials.  There  are  17  examples  in  all  when  the  four  inert  materials  are  included. 
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2.  METHODOLOGY 


Essentially,  in  Reference  1,  the  microscopic  particle  velocities  analyzed  in  Reference  2 
(Chapter  3)  were  considered  equivalent  to  the  macroscopic  particle  velocities,  Up,  observed 
experimentally  in  shocked  solid  materials  [5].  Consequently,  the  critical  particle  velocity,  Vf, 
for  lattice  break-up  and  disintegration,  derived  in  Chapter  3  of  Reference  2,  could  be  related  to  a 
corresponding  shock  v^'ave  velocity,  Usf,  and  shock  pressure,  Psf,  by  the  following  equations: 

Usf  =  Co  +  S*Vf  (1) 

Psf=PoUsfVf  (2) 

Equation  1  simply  states  the  experimentally  observed  fact  that  the  shock  wave  velocity,  Usf, 
is  a  linear  function  of  the  particle  velocity.  Up,  for  many  shock-loaded  materials  [5,  Eq.  7],  Co 
and  S  are  both  constants  and  have  been  experimentally  determined  for  a  large  number  of 
materials  [5]. 

The  shock  velocity  may  be  a  nonlinear  function  of  the  particle  velocity  for  certain  sub¬ 
stances.  In  that  case,  Us  can  be  found  from  tables  or  graphs  of  Us  as  a  function  of  Up.  Table  1 
contains  Co  and  S  information  for  the  shocked  explosives  and  materials  considered  in  this  re¬ 
port.  The  sources  of  this  information  are  also  listed. 

Equation  2  is  the  well  known  expression  for  shock  pressure  [5,  Eq.  2].  Psf  is  the  product  of 
the  material  density  (po),  shock  velocity  (Usf),  and  the  particle  velocity  (Vf). 

When  Vi  is  k-’.cw'n,  Usf  and  Psf  may  be  computed  and  compared  with  experimental  data 
(Up,  Us,  Ps)  co,:,:sponding  to  any  phenomena  which  could  be  attributed  to  lattice  disintegration 
or  bond  breakage.  This  is  basically  the  approach  followed  in  Reference  1  and  in  the  present 
report. 

The  following  relation  ,or  k  f  was  deduced  by  Fitzgerald  [2,  Chapter  3]  from  somew’hat 
unusual  particle  dynamics  concepts.  He  called  this  the  phonon  fission  velocity,  which  is: 


Vf 


Vl  ♦  yvj  t  4  •  Vi  •  Q 
2 


(3) 


It  is  approximately  equal  to  the  geometric  mean  of  Vi  and  C*  since  Vf  is  nomally  much 
less  than  Cs  so  that: 


Vf  «  M  *  Q 


(4) 


The  velocities  Vj  and  Cs  are  defined  as  follows: 


Vi 


h 

2  m  di 


(5) 


2 


Vj  =  The  De  Broglie  particle  velocity  which  is  the  limiting  free  particle  veloc¬ 

ity  whi’wh  can  occur  without  permanent  lattice  deformation  (plastic  flow);  or  the  limit  propaga¬ 
tion  velocity  for  particle-momentum  waves  in  a  stationary  lattice.  See  References  6  and  68  for 
additional  information. 


Q 


(6) 


=  A  mean  sound  velocity  defined  such  that  the  cohesive  energy,  D,  per  atom  of 
mass,  m,  is  m  Q^.  See  Appendix  B  for  additional  information. 

The  quantities  appearing  on  the  right-hand  side  of  Equations  5  and  6  are  defined  as  fol¬ 
lows: 


h  =  Planck’s  constant 


=  6.6262  X  10-27 


(gram)  (cm^) 
sec 


di  =  Closest  distance  between  the  atoms  in  a  crystal  lattice,  or  the  atomic  spacing  in 
a  slip  direction.  Units  are  angstrom  u.nits,  A  (A  =  10~®  cm). 

m  =  Mass  of  one  atom,  grams. 


Q  =  Elastic  transverse  or  shear  wave  velocity,  cm/n  -  sec. 

Cl  =  Elastic  longitudinal  wave  velocity,  cm/p  -  sec. 

A  more  general  expression  for  Q  is  given  in  Reference  3  for  materials  which  do  not  exhibit 
uniform  accoustical  wave  propagation  characteristics  in  three  orthogonal  directions.  This  rela¬ 
tion,  from  Equation  9  of  Reference  3  is: 


Cs  = 


Q^la  ^i2b  Cl2  +  Q^2a  "*■ 


•t2b 


C2 


L3 


C.2 


t3a 


C2 


18 


(7) 


Note  that  Equation  7  reduces  to  Equation  6  if  the  three  orthogonal  longitudinal  wave  velo¬ 
cities  (Cli,  Cl2.  and  Cl3)  are  equal,  and  the  six  shear  wave  velocities  (Qu,  Qib>  Q2a,  Ct2b,  Coa. 
Ct3b)  are  equal. 

It  is  maintained  in  Reference  3  that  if  m  C*^  is  the  bond  energy  per  atom,  D,  then  (or  the 
sum  of  the  squares  of  the  six  shear  wave  velocities  and  the  three  longitudinal  wave  velocities) 
must  be  an  invariant  quantity.  The  analysis  in  Reference  3  did  not  prove  that  D  =  mCj^,  but  it 
did  show  that  Cs^  was  an  invariant,  which  is  a  necessary  condition  for  this  relation  to  be  valid. 

A  more  general  form  of  Vf  was  derived  in  Reference  2,  Equation  3.17,  which  is: 

„  V,  +  /V|  ♦  4  •  V,  •  Vu  (8) 

V,  .  - ^ - 


3 


where 


V^f  =  A  velocity  defined  such  that  the  lattice  bond  dissociation  energy  per  atom  of 
mass.  m.  is  k*  m*  where  k  is  a  constant. 

That  is,  in  certain  cases,  Vu  may  not  be  equal  to  Cj  as  defined  by  Equations  6  and  7.  One 
likely  candidate  for  Vy  could  be  Cb,  the  bulk  elastic  wave  velocity.  See  Reference  47  and  Ap¬ 
pendix  B. 

It  w'as  empirically  shown  in  Reference  2  for  18  metals,  that  the  dissociation  energy  per 
atom,  D,  was  equal  to  m  This  same  idea  (that  D  =  m  Q^)  was  considered  (without  proof) 
to  be  applicable  to  polymeric  materials  in  the  analysis  reported  in  Reference  1  and  in  the  present 
study  also.  The  favorable  comparisons  of  the  computed  and  experimental  reactive  conditions 
documented  in  Reference  1  and  the  present  report  arc  consistent  with  this  assumption.  The  cred¬ 
ibility  of  this  assumption  is  also  enhanced  via  information  found  in  Reference  69.  Reference  69 
contains  estimates  of  the  energy  per  gram  (E's  =  D/m)  required  for  complete  vaporization  of 
Comp-B  and  PBX-9404.  For  both  explosives,  compared  reasonably  well  with  Cj  (Table 
B-2). 

Fitzgerald’s  derivation  of  Vf  is  based  on  the  premise  that  the  De  Broglie  relation  ( Up  = 
h.'  (2m  d'l))  must  be  satisfied  for  particle  velocities  greater  than  Vj.  This  led  him  to  the  unusual 
concepts  of  reverse  lattice  motion  and  the  extraction/depletion  of  cohesive  energy  to  create  this 
motion.  Consequently,  this  means  that  the  De  Broglie  relation  must  be  satisfied  on  a  microscop¬ 
ic  scale  within  the  shock  front  thickness  (a  few  angstroms)  by  either  reverse  lattice  motion  or 
some  other  phenomenon  even  more  bizarre! 

This  is  because  the  one-dimensional  bulk  shock  compression  [59,  pp.  463-465]  is 
definitely  not  sufficient  to  satisfy  the  De  Broglie  relation  demands  on  the  lattice  spacing  (d'j 
=  h/(2m  Up)).  The  bulk  unidimensional  shock  compression  is: 

A  -  il  -  £o  _  Us -Up 

di  Vo  Q  Us 


so  that; 


di'  <  <  d  when  Up  >  Vj  (10) 

That  is,  the  bulk  shock  compression  from  di  to  d  is  much  less  than  from  di  to  di'.  In  fact,  the 
bulk  shock  compression  (di  to  d)  is  practically  insignificant  compared  to  the  large  change  from 
di  to  di'. 

Perhaps  the  basic  reason  a  shock  wave  (or  shock  front  thickness)  is  formed  is  because  De 
Broglie’s  relation  must  be  satisfied  for  particle  velocities  greater  than  Vi.  The  shock  front  thick¬ 
ness  does  indeed  seem  to  be  very  dense  (or  more  dense  than  material  fore  and  aft  of  the  wave 
front).  This  is  observed  from  photographs  of  shocks  in  transparent  materials.  Consequently,  it  is 
plausible  to  suggest  that  it  is  in  the  shock  &ont  thickness  where  the  De  Broglie  velocity  —  wave 
length  relation  is  satisfied  (possibly  during  times  on  the  order  of  10"^^  or  10"^^  second).  Subse¬ 
quently,  there  would  be  material  relaxation  from  di'  to  d  (bulk  compression)  and  eventually  to  di 
as  relaxation  continues. 
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The  remarks  in  the  previous  paragraph  are  believed  to  apply  not  only  to  solids,  but  to  liq¬ 
uids  and  gases  as  well.  Recent  Computational  Molecular  Dynamics  (CMD)  numerical  modeling 
results  for  shock  loaded  particles  (atoms  and/or  molecules)  are  repored  in  References  63 
through  67.  It  is  recommended  that  numerical  experimentation  with  Fitzgerald’s  LPD  Vf  con¬ 
cept  be  included  in  future  CMD  studies.  The  LPD  cohesive  energy  dissipation  is  visualized  as  a 
damping  influence  which  could  possibly  explain  viscous  effects  within  the  shock  front  thickness 
[70]  and  the  shock  variable  (Ps,  Us,  Up)  decay  observed  as  a  function  of  distance  traveled. 

The  important  relationship  of  Vi  (the  De  Broglie  particle  velocity  in  an  undisturbed  or 
immobile  lattice)  with  respect  to  the  Hugoniot  Elastic  Limit  (HEL)  particle  velocity  (Uphel) 
observed  in  Elastic-Plastic  Shock  (EPS)  experiments  is  discussed  in  Reference  6.  It  is  shown 
in  Reference  68  that  critical  particle  velocities  causing  lead  azide  detonation  are  comparable  to 
Vi,  based  on  the  nitrogen  atom  mass  and  nitrogen  bond  lengths  for  this  compound. 

This  illustrates  that  refinements  with  respect  to  the  computation  of  Vi  are  feasible.  Particu¬ 
lar  atoms  or  combinations  of  atoms  and  their  associated  minimum  slip  distances,  di,  may  war¬ 
rant  special  attention.  For  example,  Fitzgerald  in  Reference  2,  Chapter  3,  computes  a  separate 
Vi  for  each  of  the  two  atoms  (Na  and  Cl)  in  salt  (NaQ).  Information  necessary  to  make  a  judi¬ 
cious  choice  of  the  atom  or  atomic  group  in  polymeric  solids  is  required  because  of  the  large 
number  of  possible  selections  which  may  not  be  physically  realistic.  Theoretical  and  experimen¬ 
tal  attempts  to  delineate  the  atoms  or  atomic  combinations  which  initiali7ft  reactions  in  explosive 
detonation  phenomena  are  documented  respectively  in  References  7  and  8.  Other  References 
reporting  similar  explorations  could  be  cited. 

The  importance  of  specific  atomic  elements  in  shocked  solids  is  very  plausible  in  view  of 
their  very  sensitive  influence  in  vibrating  bone  structure  [Appendix  C]. 

When  the  particle  velocity,  Up,  equals  or  exceeds  Vf,  self-sustained  particle  impact  or 
fission  is  possible  [2].  Fragments  of  the  broken  lattice  will  strike  other  lattice  particles,  and 
the  process  is  repeated.  Individual  atoms  could  also  be  freed  from  their  bonds  and  in  turn  strike 
other  atoms  or  lattice  fragments.  Such  an  atom  is  called  a  “free  atom”  by  Fitzgerald  [4].  The 
travel  of  a  “free  atom”  in  the  atomic  lattice  should  satisfy  the  De  Broglie  demands  even  though 
the  bulk  particles  of  shocked  materials  may  not  do  so.  See  above  discussions  of  Equations  9  and 
10  implications. 

These  “free  atoms”  in  a  shocked  solid  are  not  generated  in  a  homogeneous  array,  but  in 
a  random  heterogeneous  manner.  In  a  reactive  or  explosive  material,  these  “free  atoms”  form 
localized  “hot  spots.”  When  the  concentration  or  number  of  hot  spots  (local  particle  impacts) 
per  unit  area  reaches  a  critical  value  (dependent  on  the  explosive),  then  the  process  is  inevers- 
ible  and  detonation  is  inevitable. 

Thus  the  LDP  can  be  the  threshold  trigger  mechanism  that  creates  heterogeneous  hot  spots 
which  in  turn  causes  detonation  in  secondary  explosives  and  observable  reactions  in  many  inert 
materials. 
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3.  APPUCATIONS 


In  order  to  perform  an  assessment  analysis  for  LDP  effects,  the  following  information  is 
required.  These  data  must  be  for  the  same  material  (i.e.,  same  composition,  same  density,  Po, 


same  manufacturing  history,  etc.). 

Item 

Purpose 

1. 

m  and  di  or  mAv  tliAV 

To  compute  Vf  via  Eq.  5 

2. 

Cl  and  Q 

To  compute  Q,  via  Eq.  6  or  7  and 

Vf  via  3 

3. 

Us  =  f(Up) 

To  compute  Usf  and  Psf  from  Vf 

4. 

Experimental  shock-induced 
reactive  information  as  a 
function  of  Up,  Us,  and  Ps 

To  compare  with  the  predicted  reactive 
conditions,  Vf,  Usf,  Psf,  and  Vi 

Even  though  a  voluminous  amount  of  shock-loaded  material  information  is  available,  it  is 
difficult  to  assemble  the  above  information  for  the  same  consistent  material.  However,  enough 
sufficiently  consistent  information  has  been  found  for  four  energetic  and  three  inert  materials  so 
that  an  LDP  criteria  analysis  could  be  made.  Table  1  lists  Us  =  f  (Up)  information  for  all  the 
materials  investigated.  If  Us  =  f  (Up)  was  not  linear,  then  the  data  (Us  vs.  Up)  was  plotted  and 
Usf  was  graphically  determined.  See  Appendix  A  for  mAV  and  diAv  results. 

Table  2  contains  m,  di,  Cl,  Q,  Vi,  Vf,  Usf,  and  Psf  information  for  the  materials  under  con¬ 
sideration.  Table  3  contains  Cl,  Q  information  for  one  material,  single  crystal  PETN. 

Table  4  lists  the  comparable  comp^  ed  predicted  and  experimental  observed  reactive  condi¬ 
tions  for  all  the  materials  under  investigation.  Tables  1  through  4  also  provide  the  sources  of  in¬ 
formation  for  the  tabulated  data. 

Figures  1  through  10  depict  the  results  of  this  investigation  and  also  indicate  information 
sources.  Reference  9  is  a  very  valuable  general  information  source  for  explosives.  References 
10  and  11  are  well  known  sources  for  the  Large  Scale  Gap  Test  (LSGT)  and  Small  Scale  Gap 
Test  (SSGT)  results,  respectively,  for  explosive  sensitivity. 

Pertinent  LDP  Vf  results  for  the  energetic  and  inert  materials  arc  discussed  in  the  following 
sections. 


A.  Energetic  Materials 
1.  Pressed  TNT 

Pressed  TNT  (po  =  1.635  g/cc)  was  considered  in  Reference  1.  Since  then,  ref¬ 
erences  13, 14,  and  15  were  discovered,  which  contained  consistent  information  for  pressed  TNT 
where  Po  =  1.648  g/cc.  Psf  (20J  kbars)  compares  very  well  with  the  experimental  reaction  re¬ 
sults  (Pr  =  24  kbars)  from  Reference  15.  Both  of  these  data  points  are  plotted  in  Figure  1 
which  also  contains  experimental  data  from  References  10, 11,  and  12.  See  Reference  1  which 
also  discusses  a  similar  figure  (Fig.  3).  References  12  and  15  utilized  flat  plate  impact  to  shock 
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the  TNT  targets,  and  this  provided  precisely  defined  load  conditions.  The  Psf  values  [1,  present 
study]  compare  very  favorably  with  the  experimental  results  of  References  12  and  15.  TTiis  rein¬ 
forces  the  conclusion,  stated  in  Reference  1,  that  the  LDP  criteria  provides  an  explanation  for  the 
observed  increase  of  the  detonation  “threshold”  shock  pressure  as  the  explosive  density  increases 
or  porosity  decreases. 

2.  PBX  -  9407 

Information  contained  in  References  9  and  16  was  sufficient  for  a  consistent 
computation  of  Vf,  Usf,  and  Psf  for  PBX-9407  (Po  =  1.608  g/cc).  From  information  in  Refer¬ 
ence  9,  a  consistent  computation  of  Vf  was  made  for  Po  =  1.78  g/cc.  However,  the  Us  =  f(Up) 
relationship  was  unknown  for  this  density.  Consequently,  to  make  an  approximate  computation 
of  Usf  and  Psf  for  Po  =  1.78  g/cc,  the  Us  =  f(Up)  relation  for  Po  =  1.608  g/cc  was  employed. 
Consequently,  it  is  believed  that  these  approximate  values  of  Us  and  Psf  for  Po  =  1.78  g/cc  are 
somewhat  smaller  than  they  should  be. 

These  Psf  results  are  shown  plotted  as  a  function  of  density  (po)  in  Figure  2.  Also  shown  in 
Figure  2  are  the  SSGT  results  from  Reference  11  and  Lindstrom’s  asymptotic  Pmin  result  from 
Reference  16. 

The  Psf  magnitudes  are  somewhat  smaller  than  the  SSGT  P^  values,  but  the  variation  with 
density  is  very  similar.  The  vertical  arrow  above  the  Psf  point  for  Po  =  1.78  g/cc  indicates  that 
this  magnitude  should  probably  be  somewhat  greater  than  shown. 

For  Po  =  1.608  g/cc,  Psf  (11.1  kbars)  and  the  Pmin  result  from  Reference  16  compare  very 
well.  Pmin  (9-35  kbars)  is  the  average  value  of  Pmin  from  the  time  to  detonation  and  distance  to 
detonation  curve  fits.  Pmin  is  the  minimum  shock  pressure  necessary  to  cause  detonation.  This 
interpretation  of  Pmin  agrees  with  Lindstrom’s  alternate  assumption  with  respect  to  Pmin  which 
is  “that  the  pressure  must  be  greater  than  some  finite  limit  for  detonation  to  occur.” 

3.  Octol 

Information  contained  in  Reference  9  was  sufficient  to  allow  a  consistent  compu¬ 
tation  of  Vf,  Usf,  and  Psf  for  cast  octol  (Po  =  1.80  g/cc).  This  value  of  Psf  (24.74  kbars)  is 
shown  plotted  in  Figure  3  where  it  compares  favorably  to  LSGT  results  for  cast  octol  [10]  and 
SSGT  results  for  pressed  octol  [11]. 

4.  PETN 

An  LDP  assessment  was  made  for  both  pressed  and  single-crystal  (S.  C.)  pen- 
taerythritol  tetranitrate  explosive  which  is  commonly  referred  to  as  PETN.  The  assessment  pro¬ 
cedure  for  the  pressed  PETN  is  the  same  as  for  the  PBX-9407,  octol,  and  pressed  TNT.  Howev¬ 
er,  for  the  S.  C.  PETN,  a  modification  was  incorporated  in  the  computations.  The  illustration  of 
this  modification  was  one  reason  for  including  S.  C.  PETN  in  this  report.  Another  reason  is  to 
show  that  LDP  offers  an  explanation  for  certain  rather  low-level  (Ps  <  40  kbars)  shock-induced 
reactions  which  have  been  experimentally  observed  in  S.  C,  PETN. 

a.  Pressed  PETN 

Information  contained  in  References  9  and  17  was  sufficient  for  a  consistent 
computation  of  Vi,  Vf,  Usf,  and  Psf  for  pressed  PETN  (po  =  1.72  g/cc),  as  shown  in  Tables  1 
and  2.  In  Figure  4,  Psf  is  graphically  compared  with  test  results  from  References  11  and  17.  Psf 
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is  somewhat  larger  than  the  SSGT  results  [11].  Psf  (21.0  kbars)  is  still  larger  than  the  exper¬ 
imental  detonation  pressure  (17.0  kbars)  from  Reference  17.  TTiis  is  considered  a  reasonably 
good  comparison.  It  should  be  noted  that  although  detonation  was  not  observed  [17]  for  Ps 
below  17.0  kbars,  perhaps  it  would  have  been  for  larger  sample  sizes. 

b.  Single-Crystal  PETN 

Shock-loading  experimentation  with  S.  C.  PETN  has  been  documented  in 
References  18  through  23.  References  18,  20,  and  21  mention  that  “anomalous”  inexplicable 
reactions  were  observed  in  S.  C.  shocked  with  40  kbars  in  the  <110>  direction.  In  fact,  Refer¬ 
ence  18  indicates  that  shock-induced  decomposition  may  have  occurred  in  two  shots  where  the 
shock  pressure  was  nominally  26  kbars.  Certain  characteristics  of  the  transducer  gauges  (x-cut 
quartz)  at  these  load  levels  prevented  a  definite  resolution  of  the  observed  anomalies.  A  strong 
reaction  at  40  kbars  was  also  apparently  observed  by  Craig  [18,  20,  21]. 

The  above  reactions  at  26  and  40  kbars  were  termed  “anomalous”  because 
they  were  inexplicable  by  conventional  reaction  rate  theory.  Apparently  it  normally  required 
about  86  kbars  to  reliably  detonate  <110>  S.  C.  PETN  [18  through  22]. 

In  order  to  clarify  this  situation,  which  had  existed  for  a  number  of  years, 
experimentation  on  S.  C.  PETO  with  22.0, 35.0,  and  42.6  kbars  of  shock  input  was  performed 
[23]  with  the  following  results: 

42.6  kbars:  Detonation  occurred.  The  sample  crystal  was  larger 
than  those  employed  by  previous  investigations.  Two  different 
shots  were  fired  and  the  quartz  gauge  pressure  records  were 
practically  identical  to  one  recorded  by  Craig  at  40.0  kbars  [23, 

Fig.  1].  Q)nsequently,  it  is  assumed  that  detonation  would  have 
been  observed  at  40.0  kbars  if  the  sample  had  been  larger. 

35.0  kbars:  Strong  light  emission  was  observed.  According  to 
Reference  23,  this  is: 

“Strong  evidence  that  we  are  observing  nonequilibrium 
chemiluminescence  due  to  shock-induced  decomposition  during 
the  initiation  process.” 

22.0  kbars:  Light  emission  was  observed,  both  for  <110>  and 
<001>  orientations.  This  is  indicative  of  some  reaction  at  this 
level  of  shock-loading. 

These  reactive  pressure  data  points  (22.0, 35.0,  and  40.0  kbars)  are 
shown  in  Figure  4  on  the  crystal  Theoretical  Maximum  Density  (TMD) 
indicator  line.  Also  shown  on  the  TMD  line  is  the  LDP  Psf  value  (25.0 
kbars)  computed  via  the  information  contained  in  Tables  1, 2,  and  3. 

The  ^Sf,  Vf)  point  is  plotted  in  Figure  5  along  with  the  experimental 
results  fiom  Reference  18. 

Note  that  because  S.  C.  PETN  is  anisotropic,  the  value  of  C*  was  computed 
via  Equation  7  as  shown  in  Table  3.  The  elastic  wave  velocities  given  by  Morris  [24]  were 
employed  in  this  computation.  Note  also  that  the  average  or  “bulk”  values  of  m  and  di  were 
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used  to  compute  Vi.  The  “bulk  Vi”  LDP  Psf  prediction  (25.0  kbars)  is  remarkably  close  to  the 
experimentally  observed  light  emission  reaction  pressure,  Pr,  level  of  22.0  kbars. 

It  was  stated  in  Section  II  and  in  Reference  1  that  the  contribution  of  individ¬ 
ual  atoms  or  atomic  groups  could  be  given  special  attention  via  the  Vi  computation.  Bond 
length  distances  between  atoms  in  the  PETN  crysta’  re  given  in  Reference  25.  However,  at 
present,  the  author  does  not  possess  sufficient  infonnation  for  a  judicious  choice  of  the  correct 
atom  or  atomic  combination  to  compute  a  more  specific  or  appropriate  Vi  than  the  “bulk  Vi” 
based  on  average  values.  The  required  information  is  being  sought,  and  when  found,  a  more 
appropriate  Vi  and  the  associated  Vf,  Usf,  and  Psf  values  will  be  computed. 

B.  Inert  Materials 

1.  PPMI 

Information  for  a  shock-loaded  polymeric  substance,  PPMI,  was  acquired  from 
References  26,  27,  and  28.  Sufficient  consistent  data  were  available  so  that  an  LDP  assessment 
could  be  performed.  In  the  particle  velocity.  Up,  region  of  interest,  the  shock  velocity,  Us,  is  a 
non-linear  function  of  Up  as  shown  in  Figure  6.  From  this  figure,  Usf  was  determined  when  Vf 
was  computed.  Psf  was  15.0  kbars  which  is  very  close  to  the  “cusp”  location  [27, 28]  of  the  Us 
=  f(Up)  relationship  that  occurs  at  about  18.0  kbars. 

Reference  27  indicates  that  shock-induced  polarization  also  begins  at  approxi¬ 
mately  18.0  kbars  for  PPMI.  Thus,  both  the  Hugoniot  cusp  and  the  initial  shock-induced  polar¬ 
ization  point  is  very  close  to  the  predicted  LDP  Psf  point. 

The  authors  of  Reference  27  also  refer  to  a  shock  Hugoniot  cusp  and  shock-in¬ 
duced  polarization  in  polymethyl  methacrylate  (PPMA),  or  plexiglas,  which  occurs  at  about  20.0 
kbars.  This  coincides  rather  closely  with  the  LDP  Psf,  predicted  in  Reference  1  for  PPMA, 
which  was  21.4  kbars. 

In  Reference  26,  Graham  presents  a  strong  case  for  a  “catastrophic  shock”  event 
which  occurs  at  a  relatively  low  shock  pressure  threshold.  Above  this  pressure  threshold,  lattice 
or  chemical  “bond  scission”  and  anomalous  electrical  phenomena  occur  in  shocked  polymeric 
solids.  Furthermore,  non-equilibrium  thermodynamic  conditions  and  process  irreversibility  ex¬ 
ist  during  and  after  the  “catastrophic  shock”  event. 

Graham  postulated  that  a  “catastrophic  shock”  event  was  occurring,  based  upon 
his  observations  of  a  large  amount  of  experimental  data.  However,  a  “trigger  mechanism”  for 
low-pressure  “catastrophic  shock”  bond  scissions  has  been  difficult  to  justify  theoretically  via 
conventional  shock  physics  theory. 

Consequently,  unconventional  concepts  such  as  Fitzgerald’s  LDP  Vf  criteria  mer¬ 
it  attention.  In  view  of  the  “predictive  power”  demonstrated  in  Reference  1  and  the  present  re¬ 
port  for  a  number  of  diverse  materials,  it  is  believed  that  Fitzgerald’s  LDP  Vf  provides  the  “trig¬ 
ger  mechanism”  for  the  low-pressure  “bond  scissions”  of  the  “catastrophic  ^ock”  event. 

2.  Iron 

Reference  29  provides  a  comprehensive  review  of  the  technical  literature  up  to 
1977  with  respect  to  shock-induced  lattice-structure  phase  changes  observed  in  solid  materials. 
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The  element,  iron,  was  the  first  material  in  which  a  shock-induced  phase  change  was  found,  and 
much  experimental  work  has  since  been  performed  to  delineate  the  phenomena  which  occurs  at 
approximately  130.0  kbars. 

Figures  7  and  8  illustrate  experimental  results  for  Us  and  Ps,  respectively,  as  a 
function  of  particle  velocity.  Up,  for  shock-loaded  iron.  Sources  [30  through  34],  for  the  exper¬ 
imental  information  are  denoted  in  these  figures. 

Iron,  under  moderate  shock-loading  below  130  kbars,  exhibits  a  two-wave  struc¬ 
ture  where  the  first  one  is  essentially  an  elastic  wave  followed  by  a  slower  plastic  deformation 
wave  front.  In  general,  this  is  called  elastic-elastic  shock  phenomena  (EPSP). 

The  elastic  wave  generated  under  EPSP  conditions  is  called  the  HEL  to  distin¬ 
guish  it  from  a  dynamic  elastic  limit  which  could  occur  from  less  severe  transient  load  condi¬ 
tions  [35].  It  is  noted  in  Reference  36,  page  191,  that  the  HEL  condition  can  produce  the  same 
effect  as  a  phase  transition.  That  is,  an  additional  wave  is  formed  [6]. 

In  order  to  perform  an  LDP  assessment  for  iron,  the  velocities  Vj  and  Vf  given  in 
Table  3  were  employed  to  compute  Usf  and  Psf.  The  double  wave  structure  requires  that  these 
computations  be  performed  via  the  following  equations  [35,  p.  256]. 

The  pressure  behind  the  HEL  wave  is: 

Phel  =  PoClUphel  (11) 

=  PoClVi  (12) 

It  was  shown  in  Reference  6  that  the  steady  state  equilibrium  value  of  the  HEL 
wave  particle  velocity,  Uphel.  was  equal  to  Vi.  The  density  behind  the  HEL  wave  is: 


Equation  15  is  a  good  approximation  because  Cl,  the  longitudinal  elastic  wave 
velocity,  is  much  larger  than  Vi  so  that  Cl/(Cl-Vi)  <•  1.0. 

The  plastic  wave  pressure  corresponding  to  the  LDP  particle  velocity,  Vf,  is: 


PSf  =  PHELUSf(Vf-UpHEL)  +  PHEL 


(16) 
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Psf  =  Phel  Usf  (Vf- Vi)  +  Phel 


(17) 


The  iron  plastic  wave  velocity  is  a  linear  function  of  the  particle  velocity  as  given 
by  Equation  1.  The  values  of  Qj  and  S,  as  recommended  by  References  33  and  34  are  Q  = 

4.63  km/sec  and  S  =  1.33.  This  equation  is  plotted  in  Figure  7.  For  Vf  =  0.22  km/sec,  Usf  = 
4.92  km/sec,  and  this  point  is  shown  in  Figure  7.  Also  shown  in  Figure  7  is  the  HEL  point 
defined  by  Cl  and  Uphel-  Note  that  these  values  of  Cl  (6.04  km/sec)  and  Uphel  (0.0140  km/ 
sec,  or  essentially  Vf,  which  is  0.0144  km/sec)  are  the  same  ones  that  Minshdl  [35]  obtained 
experimentally  for  the  thick  iron  specimen  [6,  Fig.  2].  The  conesponding  Phel  via  Equation  12 
is  6.8  kbars. 


For  po  =  7.84  g/cc,  then  Phel  =  7.86  g/cc  via  Equation  14.  All  the  information 
required  to  compute  Psf  via  Equation  17  is  now  available.  This  computation  yields  Psf  =  86.3 
kbars.  The  (Psf,  point  is  shown  in  Figure  8. 

The  predicted  LDP  initiation  Psf  location  is  also  shown  in  Figure  9,  which  depicts 
experimental  data  for  the  relative  electrical  resistance  (R/Ro)  of  shock-loaded  iron  as  a  function 
of  the  shock  pressure,  Ps.  These  experimental  data  are  from  References  37  and  38.  The  earlier 
results  (1962)  from  Reference  38  show  a  rather  abrupt  jump  in  electrical  resistance  at  150  kbars, 
which  is  somewhat  beyond  the  130  kbar  transition  point  The  later  measurements  (1968),  as 
described  by  the  authors  of  Reference  37,  show  a  sudden  increase  (-  30%)  in  resistance  at  -  80 
kbar,  which  is  absent  in  the  data  of  Fuller  and  Price  [38]. 

The  authors  of  Reference  37  suggest  and  discuss  three  possible  causes  of  this 
resistance  jump  at  approximately  80  kbars.  Briefly,  these  suggested  mechanisms  are: 

1.  The  a  -*  £  transition  starts  at  a  lower  pressure  than  130  kbars. 

2.  The  shock  creates  lattice  defects. 

3.  Part  of  the  a  -  iron  is  transformed  into  a  new  phase  under  the  influence  of 
the  very  high  shear  itresses  associated  with  the  shock.  This  explanation  was 
favored  by  the  authors  of  Reference  37. 

As  shown  in  Figure  9,  the  predicted  LDP  initiation  shock  pressure  (Psf  =  86 
kbars)  is  remarkably  close  to  the  shock  pressure  region  (70  to  80  kbars)  where  the  so-called  80 
kbar  resistance  increase  occurs.  Consequently,  it  is  strongly  suspected  that  initial  LDP  is  the 
reason  for  the  resistance  rise.  LDP,  if  indeed  it  is  occurring,  could  be  responsible  for  any  of  the 
three  causative  mechanisms,  or  any  combination  of  these  mechanisms,  which  were  mentioned 
above.  LDP  would  certainly  cause  heterogenous  lattice  defects  and  could  be  responsible  for 
partial  a  -»  e  transformations  and/or  some  other  new  transient  phase  as  well. 

LDP  by  definition  implies  lattice  bond  breaking,  or  scission,  which  is  a  lattice 
defect.  When  some  (or  all)  of  these  broken  bonds  recombine  after  the  shock  passage,  they  may 
well  combine  in  the  e  lattice  form.  As  the  shock  level  (or  particle  velocity  level)  increases,  it  is 
certainly  physically  plausible  that  bond  breaking  and  recombination  in  the  e  -  structure  would 
intensify  and  even  escalate  to  a  complete  a  c  phase  transformation. 

This  postulated  LDP  description  for  the  iron  a  -*  e  transformation  is  somewhat 
analogous  to  remarks  made  in  Reference  1  with  respect  to  shock-induced  thermocouple  re- 
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sponse  and  detonation  of  PBX-9404.  The  trend  of  the  thermocouple  response  in  shocked 
PBX-9404  and  PMMA  is  similar  to  the  trend  of  the  resistance  change  in  shocked  iron  [37] 
That  is,  the  electrical  signal  deviates  from  the  expected  values  and  escalates  as  Ps  or  TJp  is 
increased. 


Consequently,  Fitzgerald’s  LDP  is  suggested  as  a  common  triggering  or  causative 
mechanism  behind  the  observed  electrical  signal  behavior  of  these  three  diverse  materials 
(PBX-9404,  PMMA,  and  iion)  under  shock-loaded  conditions. 

It  is  also  suggested  that  the  electrical  resistance  of  shocked  iron  be  experimental¬ 
ly  evaluated  again  (40  s  Ps  s  180  kbars)  where  APs  between  data  points  is  about  5  kbars  in  the 
region  between  70  and  90  kbars.  Shocked  thermocouple  output  for  these  conditions  would  be  of 
interest  also. 

3.  NaCl 

Vf,  Usf,  and  Psf  were  computed  for  S.  C.  NaCl,  which  was  shock-loaded  in  the 
<100>  direction.  The  Us  versus  Up  relationship  employed  was  from  Reference  39  (Table  1). 
This  relation  (Us  =  3.528  +  1.343  *  Up  km/sec)  was  a  fit  to  Los  Alamos  Scientific  Laboratory 
(LASL)  data  obtained  at  “low”  pressures.  A  plot  of  this  data  in  Figure  4  of  Reference  39  reveals 
that  <100>  -  cut  S.  C.  samples  were  employed. 

The  Cl  and  Q  values,  to  compute  Cs,  were  from  Reference  40.  An  average  value 
of  m  and  di  was  employed  to  compute  Vj.  This  value  of  di  (2.82  A )  agrees  rather  closely  with 
similar  information  given  in  textbook  examples  [41,  pp.  652-655]  for  the  NaCl  crystal  lattice 
structure.  An  average  value  of  m(mAv)  is  considered  appropriate,  since  shock-loading  it  in  the 
<100>  direction  would  involve  both  Na  and  Cl  atoms  whose  masses  are  not  greatly  different 
(Appendix  A). 


Utilizing  the  information  described  above  and  given  in  Tables  1  and  2,  Vf  and  Psf 
for  <100>  cut  S.  C.  NaCl  was  computed  to  be: 

Vf  =  0.253  km/sec 

Psf  =  21.0  kbars 

Reference  29  discusses  the  low-pressure  (20  -  30  kbars)  phase  transition  contro¬ 
versy  with  regard  to  NaQ.  Royce  [40]  states  that  a  “weak”  transformation  is  observed  at  Pst  = 
29  kbars  in  <100>  -  cut  S.  C.  NaCl.  TTie  two-wave  velocities  differed  by  only  1/2  to  1%,  and 
rather  sensitive  quartz  transducers  must  be  employed  to  detect  it.  In  Reference  42,  Weidner  and 
Royce  conclude  that  “the  29  kbar  discontinuity  for  100  crystals  can  be  considered  a  second  HEL 
caused  by  strain  hardening  on  the  major  slip  system.”  The  above  Psf  value  (21  kbars)  is  indica¬ 
tive  that  LDP  could  begin  a  few  kilobars  below  the  observed  Pst  point,  and,  as  such,  also  be  the 
responsible  mechanism. 

Indicator  lines  for  Psf  and  Pst  are  drawn  in  Figure  10  where  shock-induced  polar¬ 
ization  maximum  amplitude  signals  (volts)  are  plotted  as  a  function  of  Ps.  These  test  results  are 
from  Reference  43,  Table  1.  The  specimens  were  <100>  oriented  single-crystals  whose  diame¬ 
ter  was  25  mm.  The  specimen  thickness,  which  influences  the  signal,  varied  as  tabulated  in  Fig¬ 
ure  10.  This  thickness  variation  is  cot  enough  to  influence  the  trend  of  the  signal  output  as  a 
function  of  shock  pressure.  In  particular,  note  the  six-fold  increase  in  signal  amplitude  from  Ps 
=  16  kbars  to  Ps  =  25  kbars.  These  specimens  had  essentially  identical  thicknesses,  so  there  is 
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definitely  not  a  thickness  effect  on  the  signal  magnitudes.  Obviously,  the  LDP  Psf  and  its  prob¬ 
able  consequence,  Pst,  are  remarkably  close  to  the  onset  and  step>-like  increase  of  the  electrical 
polarization  test  signals. 


4.  DISCUSSION 

In  Reference  1  and  the  present  study,  a  total  of  nine  different  explosives  and  four  different 
inert  materials  have  been  investigated  for  possible  LDP  effects.  This  has  been  accomplished  by 
computing  Vf  and  Psf  for  these  materials,  and  comparing  Psf  with  the  shock  pressures  where 
some  type  of  shock-induced  reaction  is  occurring.  On  the  basis  of  these  favorable  comparisons 
it  is  suggested  that  the  following  shock-loaded  material  phenomena  is  explained  by  Fitzgerald’s 
LDP-critical  particle  velocity,  Vf,  concept  as  applied  to  shocked  solid  substances: 

1.  Threshold  detonation  pressure,  Pd,  for  explosives.  That  is,  Pd  -  Psf.  [1,  this  report]. 

2.  Threshold  detonation  pressure,  Pd,  variation  with  density  (po).  [1,  this  report]. 

3.  Roth’s  observation  [45]  that  although  Pd  for  explosives  varied  with  density  (po),  the 
particle  velocity  where  Pd  occurred  was  essentially  invariant  (Table  5,  Fig.  11). 

4.  Experimentally  observed  “heterogenous  hot  spots”  in  shocked  explosives.  [2,  this 
report]  The  bond  or  lattice  break-up  is  not  always  homogeneous,  and  will  not  neces¬ 
sarily  occur  all  at  once  at  the  same  time. 

5.  A  “catastrophic  shock”  event  occurring  at  low  shock  pressures  where  “bond  scission” 
occurred.  Such  an  event  was  suggested  by  Graham  [26]  to  account  for  effects  which 
have  been  observed  from  examining  a  large  amount  of  experimental  data.  These  effects 
(and  the  LDP  assessment  examples)  are: 

a.  Anomalous  electrical  behavior  [1,  this  report]. 

b.  Increased  chemical  activity  [this  report]. 

c.  Polymorphic  phase  changes  [this  report]. 

In  the  discussion  of  item  3,  Table  5  and  Figure  11  were  mentioned  to  illustrate  that  Vf  is 
practically  invariant  for  these  energetic  materials.  Table  6  and  Figure  12  were  prepared  for  the 
inert  materials  considered  in  Reference  1  and  the  present  study.  Vf  range  indicators  are  also 
shown  for  the  metallic  element  data  given  in  Reference  2,  Chapter  3,  Table  3.1.  All  of  these  Vf 
values  are  somewhat  less  than  one  km/sec,  which  is  indicative  of  relatively  low  Psf  values,  com¬ 
mensurate  with  item  5  above. 

An  important  spin-off  from  the  LDP  investigations  has  been  the  results  reported  in  Refer¬ 
ence  6.  This  document  clearly  shows  the  important  relationship  between  Vf  (the  De  Broglie 
particle  velocity)  and  classic  experimental  IffiL  particle  velocity,  UpHEL-  Since  the  De  Broglie 
velocity  is  important  at  EPSP  conditions,  it  is  reasonable  to  assume  that  it  applies  to  condi¬ 
tions  where  the  particle  velocity  is  greater  than  Vi  or  Uphfi  .  In  fact,  Fitzgerald’s  LDP  Vf  deri¬ 
vation  is  based  on  the  premise  that  the  De  Broglie  relation.  Up  =  h/(2mdi'),  must  be  satisfied 
for  Vj  i  Up  s  Vf. 
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The  De  Broglie  relation  must  also  be  satisfied  for  Up  »  Vf,  where  no  more  reverse  motion 
of  the  lattice  can  occur,  since  its  lattice  bonding  energy  is  depleted.  Consequently,  something 
has  to  occur  (a  polymorphic  phase  change,  atomic  spacing  compression,  melting,  vaporization, 
or  some  other  phenomena)  so  that  Up  =  h/(2mdi')  is  satisfied.  The  stiffened  Hugoniot  relation¬ 
ship  (Us  vs.  Up)  observed  at  the  higher  shock  pressures  (or  high  Up)  for  certain  materials  is  per¬ 
haps  indicative  of  the  severe  atomic  spacing  compression  (di')  requirements,  which  may  apply 
only  to  “  free  atoms”  and  not  to  the  bulk  particles. 

Quite  frankly,  the  author  is  amazed  at  the  “predictive  power”  of  Fitzgerald’s  LDP  Vf  con¬ 
cept  as  applied  in  Reference  1  and  the  present  investigation  to  shock-loaded  solid  materials.  It 
would  be  rather  difficult  to  devise  a  strictly  empirical  procedure  to  predict  as  well  for  so  many 
diverse  substances.  In  direct  contrast  to  an  entirely  empirical  procedure,  the  LDP  Vf  concept 
has  a  considerable  amount  of  unconventional  rationale  behind  it  (reverse  lattice  motion  via  lat¬ 
tice  bond  cohesive  energy  utilization).  The  agreement  between  the  LDP  Vf  predictions  and  the 
experimental  reactive  conditions  for  seventeen  different  cases  can  hardly  be  considered  fortu¬ 
itous. 

Consequently,  a  “bottom  line”  result  of  this  investigation  is: 

The  LDP  Vf  effect  is  strongly  believed  to  provide  a  “common  cause”  or  “trigger 
mechanism”  for  the  relatively  low-pressure  shock-induced  reactive  behavior  ob¬ 
served  experimentally  in  certain  solid  materials. 


5.  RECOMMENDATIONS 

Some  recommendations  for  additional  work  have  already  been  made,  either  explicitly  or 
implicitly,  in  the  previous  sections.  For  easy  reference,  these  are  restated  here. 

A.  The  importance  of  the  De  Broglie  momentum-wave  length  relation  with  respect  to 
particle  velocity  propagation  in  shock-loaded  solid  materials  is  apparent  from  References  1,  6, 
68,  and  the  present  report.  The  De  Broglie  relation  implications  should  be  addressed  in  future 
shock  physics  research,  particularly  via  numerical  CMD  experimentation. 

B.  Additional  available  experimental  shocked  material  data  should  be  investigated  for 
possible  LDP  Vf  effects.  Certain  anomalous  behavior  observed  at  relatively  low-level  shock¬ 
loading  could  possibly  be  explained  as  a  Vf  effect,  or  as  a  Vi  effect  [6, 68]. 

C.  The  electrical  resistance  of  shocked  iron  should  be  experimentally  evaluated  again 
(40  <  Psf  <  180  kbars)  where  the  APs  between  data  points  is  no  greater  than  5  kbars  in  the  region 
between  70  and  100  kbars.  Shocked  thermocouple  output  for  these  conditions  would  also  be 
desirable. 

D.  Additional  experimentation  (polarization,  thermocouples,  electrical  resistance,  etc.) 
should  be  performed  on  NaCl  where  Ps  is  varied  from  about  15  to  35  kbars  in  APs  increments 
no  greater  than  5  kbars. 
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TAELE  1.  Shocked  Material  Information. 


Material 

Po 

grams/cm^ 

Co 

cm/p-sec 

S 

Source 

PBX-9407 

1.60 

(pressed) 

0.1328 

1.993 

[16] 

Octol 

1.803 

(cast) 

0.221 

2.51 

[9] 

TNT 

1.648 

(pressed) 

0.238 

2.215 

[14] 

averaged  results 

PETN 

1.72 

(pressed) 

0.183 

3.45 

[17] 

1.778 

(crystal) 

0.274 

1.81 

[22] 

Iron 

7.84 

0.463 

1.33 

[33,  34] 

NaCl 

2.165 

0.3528 

LASL 

Data  Fit 

1.343 

LASL 

Data  Fit 

[39] 

PPMI 

1.425 

Non-linear  11$  versus  Up 
variation.  (Fig.  6) 

[26,  28] 
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TABLE  2.  Compilation  of  Vj,  Q,  Vf,  Usf,  and  Psf  Results. 
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Usf  for  PBX-9407  (po  =  1.78  G/CC)  was  computed  from  data  for  Po  =  1.608  G/CC. 


TABLE  3.  Single  Crystal  PEW  Elastic  Wave  Velocities  [24], 


X  Type  Mode  Direction  of 


Motion 

—  A 

km/sec 

1 

L 

100 

3.1140 

2 

T 

001 

1.6841 

3 

T 

010 

1.4893 

4 

L 

001 

2.6167 

5 

T 

010 

1.6822 

6 

T 

100 

1.6822 

•7 

/ 

L 

110 

2.9308 

8 

T' 

1 

-no 

1.8197 

9 

T 

001 

1.6866 

=  42,0034 

(km/sec)" 

= 

Scj  _ 
18 

2.3335  (km/sec)" 

Cs  =  1.5276  km/sec  via  Eq.  7 
L  =  Longitudinal  wave,  Cl 
T  =  Transverse  or  shear  wave,  Q 


C  ‘ 

(km/sec)- 

9.6970 

2.8362 

2.2180 

6.8471 

2.8298 

2.8298 

8.5896 

3.3113 

2.8541 
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TABLE  4.  Theoretical  and  Experimental  Reactive  Shock  Information  Tabulation. 


Material 

Usf 

Psf 

Exp.  Data 

Exp.  Data  Remarks  with  Regard  to 

cm/p-sec 

kbars 

kbars  or 

Experimental  Data 

cm/p-sec 

Item 

Ref. 

PBX-9407 

0.008327 

11.10 

9.35 

Pmin 

[16] 

Average  result  from  time  and 

po=  1.60 

0.03430 

distance  to  detonation  test  data. 

g/cc 

0.2012 

See  main  text. 

17.00 

Pso 

[11] 

SSGT  50%  firing  stimulus. 

(Fig-  2) 

Octol 

0.008764 

24.74 

25.00 

P50 

[11] 

SSGT  50%  firing  stimulus 

o 

00 

(1 

o 

Q. 

0.04208 

of  pressed  octol.  (Fig.  3). 

g'cc  (cast) 

0  3266 

22.00 

Pi 

[10] 

LSGT  detonation  threshold  Ps. 

Cast  75.25  octol. 

Po  =  1.79  g/cc 

TNT 

0.008321 

20.49 

24.00 

Pr 

[15] 

Reaction  noted  in  impact  test  of 

Po  =  1.648 

0.03847 

pressed  TNT. 

g'cc  (pressed)  0.3232 

Po  =  1 .648  g/cc 

PETN 

0.008352 

21.07 

17.0 

Pd 

[17] 

Detonation  observed.  No 

Po=  1.72 

0.03903 

detonation  at  15.0  kbars  in 

g/cc  (pressed)  0.3176 

sample. 

PETN 

0.008437 

24.87 

42.6 

Pd 

[23] 

Detonation  observed. 

Po=  1.778 

0.04037 

g'cc  (crystal) 

0.3470 

35.0 

Pr 

[23] 

High  level  reaction  observed. 

22.0 

Pr 

[23] 

Some  reaction  observed. 

PPMI 

0.009032 

15.13 

18.0 

Pr 

[27] 

Cusp  in  Us  vs.  Up  relation, 

Po=  1.425 

0.03677 

plus  start  of  shock-induced 

g/cc 

0.2870 

polarization. 

Iron 

0.00144 

86.26 

70.0 

*v 

[35] 

Increase  in  electric  resistance 

Po  =  7.84 

0.02197 

to  80.0 

measurement.  (Fig.  9). 

g/cc 

0.4923 

130.0 

Pr 

[30] 

Phase  transition  observed. 

NaCl 

0.002421 

21.16 

16.0 

Pr 

[43] 

Shock-induced  electrical  signal 

p0  =  2.165 

0.0253 

to  25.0 

increase.  (Fig.  10). 

g/cc 

0.3867 

<100>  cut 

single  crystal 

29.0 

Ps. 

[40] 

Phase  transition  observed. 
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TABLE  5.  Summary  of  LDP  Assessment  Results  for  Explosives. 


Explosive 

Po  (g/cc) 

Vi  (km/sec) 

Vf  (km/sec) 

Psf  (kbars) 

Octol 

1.80 

0.08764 

0.4208 

24.74 

Comp-B3 

1.70 

0.08571 

0.4170 

26.59 

PBX-9407 

1.78 

0.08614 

0.4148 

15.91 

PBX-9404 

1.842 

0.08879 

0.4109 

26.16 

Comp)-B 

1.70 

0.08228 

0.4077 

25.16 

PETN  (S.C.) 

1.778 

0.08437 

0.4037 

24.87 

H-6 

1.75 

0.07959 

0.3903 

23.87 

PETN  (P) 

1.72 

0.08352 

0.3903 

21.07 

TNT(P) 

1.648 

0.08321 

0.3847 

20.49 

TNT(P) 

1.635 

0.08299 

0.3739 

18.09 

PBX-9407 

1.608 

0.08327 

0.3430 

11.10 

Tetryl 

1.70 

0.07729 

0.3411 

17.00 

Tetryi 

1.50 

0.07443 

0.3340* 

8.53* 

ViAV  =  0.08302  km/sec  •  Correction  of  Vf,  Psf  given  in  [1] 

VfAV  =  0.3871  km/sec 

(P)  =  Pressed 

(S.C.)  =  Single  Crystal 
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TABLE  6.  Summary  of  LDP  Results  for  Inert  Materials. 


Substance 

Po  (g/cc) 

Vi  (km/sec) 

Vf  (km/sec) 

Psf  (kbars) 

PMMA 

1.18 

0.1417 

0.5157 

21.0 

PPMI 

1.425 

0.09032 

0.3677 

15.0 

NaCl  (S.  C.) 

2.165 

0.02421 

0.253 

21.0 

Iron 

7.84 

0.01441 

0.2197 

86.0 

20  Metals 

0.072 

[2] 

to 
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COMPUTATION  OF  Mav  AND  Diav  FOR  CHEMICAL 
COMPOUNDS  AND  MIXTURES  OF  COMPOUNDS 


Certain  solid  materials  considered  in  this  study  were  either  chemical  compounds  (PETN, 
PPMI,  and  NaCl),  or  mixtures  of  chemical  compounds.  For  these  compounds  and/or  mixtures, 
the  weighted  average  mass,  mav>  of  a  single  atom  in  the  material  was  desired. 

First,  it  was  necessary  to  compute  the  mass  of  a  single  atom  for  each  of  the  elements  con¬ 
tained  in  the  solid.  Each  solid  was  composed  of  one  or  more,  but  not  all,  of  the  following  ele¬ 
ments: 

Carbon,  C;  Hydrogen,  H;  Nitrogen,  N; 

Oxygen,  O;  Sodium,  Na; 

Chlorine,  Cl;  Flourine,  F. 

The  mass  of  a  single  atom  of  these  elements  is: 


MW  Grams/ Gram  -  Mole  Grams 

Nav  Atoms/ Gram  -  Mole  Atom  (A-1) 

Where: 


MW 

Nav 


Molecular  Weight, 


Grams 

Gram  -  Mole 


Avogadros  Number 


6.02252  X  10*23 


Atoms 
Gram  -  Mole 


Table  A-1  lists  MW  and  m  for  each  of  the  elements  in  the  above  list.  Values  of  Nav  and 
MW  are  from  various  chemistry  text  books  and  handbooks. 

To  compute  the  average  weight  (mav)  of  an  atom  in  the  material,  the  chemical  formula  or 
proportional  chemical  composition  is  required.  Of  course,  the  weight  (m)  of  each  elemental 
atom  must  also  be  known,  since  miav  is  just  a  weighted  average  of  the  elemental  atoms  in  the  ma¬ 
terial. 


When  mav  is  computed,  then  the  average  space  between  atoms  (diav)  is  given  by  the  follow¬ 
ing  relation: 


diav 


(A-2) 


Computations  of  mav  and  diav  for  PBX-9407,  octol,  PETN,  PPMI,  and  NaCl  are  included 
in  this  appendix. 
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TABLE 

A-1.  Mass  of  a  Single 

Element 

MW 

- 

Grams 

Gram-Mole 

Carbon,  C 

12.011 

Hydrogen,  H 

1.008 

Nitrogen,  N 

14.008 

Oxygen,  0 

16.00 

Chlorine,  Cl 

35.45 

Flourine,  F 

18.99 

Sodium,  Na 

28.99 

Atom  for  Selected  Elements. 


Nav 

m 

Atoms 

Grams 

Gram-Mole 

Atom 

6.02252(1023) 

1.9943(10-23) 

6.02252(1023) 

0.1674(10-23) 

6.02252(1023) 

2.3259(10-23) 

6.02252(1023) 

2.6567(10-23) 

6.02252(1023) 

5.8874(10-23) 

6.02252(1023) 

3.1546(10-23) 

6.02252(1023) 

3.8173(10-23) 

£BX=9.4Q7: 


po  =  1.608  grams/cm^ 


Chemical  composition:  C1.41  H2.66  N2.54  O2.54  CI0.07  F0.09 

C1.41:  1.41  •  1.9943(10-23)  =  2.8120(10-23)  grams 

H2.66:  2.66  *  0.1674(10-23)  =  0.4453(10-23)  grams 

N2.54:  2.54  *  2.3259(10-23)  =  5.9078(10-23)  grams 

O2.54:  2.54  *  2.6567(10-23)  =  6.7480(10-23)  grams 

CI0.07:  0.07  *  5.8874(10-23)  =  0.4121(10-23)  grams 


F0.09: 


9.31  atoms  16.6091(10-23)  grams 


.  16,62^  ^  l,7840(l(r-)£2!2i 

9.31  atom 

at,,  =  .  „.0,45(ir«)cm^ 

l.OUo 

dlAV 


2.2303(10-«)cm 


Octol 


Po  =  1.80  grams/cc 

Chemical  comnnsitinri-  C,  •70  TTn  CO  25  O''  C5%  HMX,  25%  TNT) 

Ci,78:  1.78  *  1.9943(10-23)  =  3.5499(10-23)  grams 

H2.58;  2.58  *  0.1674(10-23)  =  0.4319(10-23)  grams 

N2.36:  2.36  *  2.3259(10-23)  =  5.4891(10-23)  grams 


O2.69: 


-23\  ^7  1. 


9.41  atoms 


16.6174(10-23)  grams 


=  .  1.7659(10-“)^ 

9.41  atom 


mAv  _  17.o592(10-^'*) 

Qo  1.80 

=  2.1407(10-®)cm 


=  9.8107(l(r24)cm 


Pentaeiythritol  Tetranitiate  (PETN) 


Po  =  1.773  grams/cc 

Chemical  composition:  C5  Hs  N4  O12 

C5:  5  *  1.9943(10-23)  =  9.9715(10-23)  grams 

Hg:  8  *  0.1674(10-23)  =  1.3392(10-23)  grams 

N4;  4*  2.3259(10-23)  =  9.3036(10-23)  grams 


-23'!  ^ 


29  atoms 


52.4947(10-23)  grams 


,  52.4947(1(^23)  ^  ,  8,o2(l(r23)  ^ 

29  ^  atom 


18.102(10-2'*)  grams 


1.773 


/warns) 


=  10.2098  *  (10-2") 


llAV  = 


2.1694(10-«)  cm 


Polypyromellitimide  (PPMI) 


po  ■  1.425  grams/cc 

Chemical  composition  [26]:  (€22  Hio  N2  05)n 
C22-.  22.0  *  1.9943  (KT^^)  =  43.8746  (KT^^)  grams 

Hio;  10.0*0.1674(10-23)  =  1.6740  (10-23)  grams 

N2:  2.0  *  2.3259(10-23)  =  4.6518  (10-23)  grams 

Qsl _ 5.0  *  2.6567  fl0-23l  =  13.2835  (10-231  grams 

39.0  atoms  63.4839  (10-23) 


mAv 


63.4839  (10-23) 
39!o 


1.6278  (10-23) 

'  '  atom 


^lAV 

dlAV 


16.278  (10-2^)  .  . 

— ^  i =  11.4231  (10-2^)  cm3 

Qo  1.425  V  ) 

2.2521  (10  -®)  cm 


Sodium  Chloride 


Po  - 

2.165  grams/cc 

Chemical  Composition:  NsCl 

Na: 

1.0  *  3.8173  (10-23)  = 

3.8173  (10-23)  grams 

CL 

1.0  *  5.8874  n0-23l  = 

5.8874  ('10-23)  pj-ams 

2.0  atoms 

9.7047  (10-2^)  grams 

9.7047  (10-23) 

(l(r23)  grams 

mAv 

=  - - -  =  4-8525 

2.0 

dlAV 

48.525 

=  22.4127  1 

Qo  2.165 

dlAV 

=  2.8195  (l0-»)  cm  -  2.82 

(lO"*)  cm 

This  value  for  di  agrees  well  with  the  atomic  separation  distance  for  NaCl  given  in  Refer¬ 
ence  41,  pages  652  through  655,  as  a  textbook  example. 
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APPENDIX  B 


RELATIONSHIPS  BETWEEN  ELASTIC  WAVE 
VELOCITIES  AND  THE  COHESIVE  ENERGY 
OF  SOUD  MATERIALS 

Since  the  publication  of  Reference  1  and  in  the  latter  part  of  this  investigation,  some  impor¬ 
tant  technical  literature  [46  through  49]  was  discovered  which  relates  the  shock  wave  and  par¬ 
ticle  velocity  parameters  (Cq  and  S  in  Eq.  1)  to  the  material  cohesive  energy.  This  relationship, 
which  was  suggested  by  McQueen  and  Marsh  [46],  has  been  thoroughly  investigated  by  Rodean 
[47, 48,  49].  The  relation,  employed  by  Rodean  in  his  extensive  comparisons  of  experimental 
data,  is: 


or 


Ci  =  2S2  Ec 

where 

Co  *«  Q,  =  bulk  sound  velocity  =  Qo 

K  =  bulk  modulus,  gram  (cm/sec)^ 

Ec  =  material  cohesive  energy  per  unit  mass,  ergs/gram  or  (cm/sec)^. 

-  Hs,  sublimation  energy  per  unit  mass,  ergs/gram  or  (cm/sec)^. 

Note  that: 

E  -  „  Uc  _  D 

‘  AW  MW  in  (B_2) 

where: 

Uc  =  cohesive  energy  per  molecular  weight  (MW,  grams)  of  a  compound  or  atomic 

weight  (AW,  grams)  of  an  element 

=  Nav  D  for  elements 

=  Ca  Nav  D  for  compounds  or  mixtures 

D  =  cohesive  energy  per  atom,  gram  (cm/sec)^ 

m  =  mass  of  an  atom  for  elements  or  average  mass  of  an  atom  for  compounds  and 

mixtures,  grams 

AW  =  Navto  =  atomic  weight,  grams 
MW  =  Ca  Nav  tn  =  molecular  weight,  grams 


Nav  =  Avogadro’s  number  [Appendix  A] 

Ca  =  constant  for  a  given  compound  or  element 

Since  S  -  1.4  -  v'2  for  many  materials,  then  Equation  b-1  becomes  approximately 


Q, 


Ec  = 


£ 

m 


or 

2/eI  =  ijD/m 


(B-3^ 


Rodean  provides  extensive  comparative  results  for  Equation  B-1  (56  metals,  12  com¬ 
pounds,  and  4  polymers). 

The  general  concept  that  the  cohesive  energy  of  a  solid  is  related  to  the  bulk  sound  wave 
velocity  is  not  new,  dating  prior  to  World  War  I.  Rodean  [47]  cites  recent  developments  con¬ 
tained  in  References  50  and  51.  Otpushchennikov  [51]  develops  the  following  relation  which  is 
also  derived  in  a  number  of  solid  state  physics  textbooks. 


2  ^  mi  ni  _  mi  ni  _D  (B-4) 

9  9  m 

where  mi  and  ni  are  exponents  in  the  potential  energy  relation.  Otpushchennikov  employed 
mi  •  ni  =  45  to  make  his  comparisons  with  experimental  data  (Cb^  =  5  Ec). 

Rodean  [49]  notes  that  mi  ni  =  32  would  provide  a  better  fit  to  the  data.  If  mi  ni  =  36, 
then  Equation  B-4  is  equivalent  to  Equation  B-3. 

In  Reference  49,  Rodean  discusses  the  results  of  Nemilov  [52]  who  derived  the  following 
relation  between  Q,  the  shear  wave  velocity,  and  the  cohesive  energy,  Ec. 


C? 


(B-5) 


Nemilov  presented  comparative  results  for  17  elemental  metals  and  two  compounds 
which  agreed  very  well  with  Equation  B-5.  He  also  presented  data  for  certain  elements  and 
compounds  where: 


C? 


2  Ec 


2  D 
m 


(B-^) 


For  Si02  (glass)  and  alkali  halide  compounds,  Nemilov’s  comparative  results  indicated 

that: 


(B-7) 


B-2 


However  for  the  alkali  halides,  Nemilov  actually  employed  Eb  (binding  energy  required  to 
transform  the  solid  to  an  ionized  gas)  instead  of  Eg  (cohesive  energy,  energy  required  to  trans¬ 
form  the  solid  to  an  un-ionized  or  neutral  gas).  J  Ec  is  employed  for  these  alkali  halides,  then 
Equation  B-5  is  satisfied  rather  well  as  shown  by  the  comparative  results  in  Table  B-1.  The  val¬ 
ues  of  Uc  for  the  alkali  halides  in  Table  B1  are  from  Reference  53  where  they  are  called  the  dis¬ 
sociation  energy.  In  these  cases,  Eb  -  2  Ec. 

Grimvald  [54,  pg.  1]  points  out  that  the  binding  energy,  Eb,  may  be  different  from  the  cohe¬ 
sive  energy,  Ec,  particularly  for  ionic  solids.  He  also  cautions  that  the  distinction  between  the 
binding  energy  and  cohesive  energy  is  not  always  made  in  the  literature.  Rodean  [49]  comments 
on  the  important  distinction  between  Ec  and  Eb- 

It  is  noteworthy  that  Fitzgerald  [55]  in  1964  had  also  suggested  Equation  B-5  and  pres¬ 
ented  comparative  results  for  12  metals.  That  is,  in  Reference  55,  C*  was  the  shear  wave  veloc¬ 
ity  and  was  favorably  compared  to  /o/mfor  metals.  Later,  by  1966  [2],  Fitzgerald  had  derived 
the  relation  for  Cj  utilized  in  the  main  body  of  this  report  (Eq.  6)  which  is: 

q2  ^  cj  +  2 

6 

This  is  essentially  a  weighted  average  of  the  sum  of  the  squares  of  the  wave  velocities  for 
isotropic  materials.  A  somewhat  more  general  result  for  certain  anisotropic  materials  is  given  by 
Equation  7  [3,  for  additional  comments  on  the  implications  of  the  general  Cs  definition]. 

For  many  materials,  Q  as  defined  by  Equation  B-8  is  practically  equal  to  Q  [1,  Table  2; 
this  report.  Table  2].  This  explains  the  good  agreement  with  /D/mfor  both  Cs  and  Q,  where 
data  for  the  same  elemental  metals  were  compared. 

It  is  also  interesting  that  both  Q  and  Ct  are  approximately  equal  to  Cp,  the  average  sound 
velocity  employed  in  Debye’s  theory  of  specific  heats.  That  is: 

Cp  •«  Ct  “  Cs  (®~9) 

Grimvald  [54]  provides  the  following  relations  for  Cp: 


(B-10) 

(B-11) 


where  v  is  Poisson’s  ratio,  which  is  approximately  0.31  for  many  solid  materials. 


Cd  is  compared  to  yD/rain  Table  B-1  for  16  alkali  halide  compounds  where  Co  was 
computed  via  Equation  B-11.  Values  of  Q  and  Q  also  are  shown  in  Table  B-1  for  some  of 
these  materials.  Similar  comparative  information  for  polymethyl  methacrylate  (PMMA)  and 
iron  is  also  presented  in  Table  B-1.  The  value  of  Uc  for  PMMA  is  from  Reference  56  by  Zhur- 
kov.  Zhurkov’s  results  are  mentioned  by  Krausz  and  Eyring  [57]  in  their  discussion  of  the  rate 
theory  of  fracture  (Chapter  4,  Section  4.4). 

Grimvald  [54]  notes  that  Co  can  be  related  to  the  Debye  temperature,  6.  He  and  Sjbdin 
[58]  provide  a  correlation  of  cohesive  energy  and  Debye  temperature  data  for  a  variety  of 
materials. 

The  information  from  different  sources  cited  in  this  Appendix  clearly  indicates  that  jD/m 
is  related  to  the  velocity  of  sound  or  elastic  wave  propagation  (Cb,  Q,  C*,  and  Co)  in  materials. 
Q,  Q,  and  Co  are  approximately  equal  in  magnitude.  Cb  is  usually  somewhat  larger  than  Q,  Cs, 
and  Cd- 

In  the  present  report  and  in  Reference  1,  Cs  was  employed  in  the  LDP  analysis  as  a  measure 
of  y  D/m  for  substances  (primarily  energetic  materials)  where  the  cohesive  energy  was  some¬ 
what  uncertain,  particularly  under  shock-loading  circumstances.  If  Q  or  Cd  had  been  employed 
for  Vu  in  Equation  8,  the  Vf  results  would  be  essentially  the  same. 

If  Cb  had  been  utilized  for  Vlt  in  Equation  8  to  compute  Vf,  then  Vf,  Usf  and  Pjf  would  gen¬ 
erally  be  larger  than  the  results  given  herein  and  in  Reference  1.  Note  that  Vf  is  approximately 
equal  to  the  geometric  mean  value  of  Vj  and  Vlt. 


Vf  ««  yvi  * 

«  yvi  *  yi; 


(B-12) 


Thus  Vf  is  not  extremely  sensitive  to  small  changes  in  Vlt  or  Eg. 

From  the  comparative  results  presented  herein  (Section  III,  A)  and  in  Reference  2,  it  is  be¬ 
lieved  that  Cs^  (or  even  or  Cd^  is  a  good  measure  of  the  cohesive  energy  for  energetic  mate¬ 
rials.  This  belief  is  strengthened  by  the  estimates  of  the  energy  per  gram  (E*  =  D/m)  necessary 
for  complete  vaporization  of  Comp-B  and  PBX-9404  which  was  found  in  Reference  69.  This 
information  is  listed  in  Table  B-2  along  with  Ct,  Cd,  C*  for  these  explosives.  Note  that 
yE^  compares  reasonably  well  with  these  elastic  wave  velocities  for  both  explosives. 

From  the  information  given  in  Table  B-3  and  plotted  in  Figure  B-1,  it  appears  that  is  a 
good  measure  of  the  heat  of  detonation  (AHdet)  of  energetic  materials.  That  is,  this  information 
indicates  that: 


eg 


A  II  _  ^  *^det 

a  xIdet  =  - 

m^y 


(B-13) 


where: 

AHdet  =  Heat  of  Detonation,  Mega  Joules/Kilogram  or  (km/sec) 
The  values  of  Cb  and  AHdet  in  Table  B-2  are  from  Reference  9. 
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TABLE  B-1.  Comparison  of  Elastic  Wave  Velocities  and  Cohesive 
Energy  for  Selected  Inert  Materials. 


Compound  MW 

Uc 

Uc/MW 
or  D/m 

VUc/MW 
or  Vb/m 

Ct 

Cd 

Cs 

grams 

gram 
(cm/sec)^ 
*  10-10 

(cm/sec)^ 
*  10-10 

km/sec 

km/sec 

km/sec 

km/sec 

CsBr 

212.8 

381.0 

1.790 

1.330 

1.16 

1.30 

KBr 

119.0 

382.8 

3.217 

1.793 

1.695 

1.90 

KCl 

74.6 

427.0 

5.7239 

2.392 

1.77 

[44] 

1.99 

2.09 

KF 

58.1 

494.0 

8.5026 

2.916 

2.59 

2.90 

KI 

166.0 

331.0 

1.9939 

1.412 

1.47 

1.65 

LiBr 

86.9 

423.0 

4.8677 

2.206 

2.07 

2.32 

Li  Cl 

42.4 

469.0 

11.0613 

3.326 

3.06 

3.43 

LiF 

25.9 

577.0 

22.2780 

4.720 

4.88 

[44] 

5.46 

3.87 

Li  I 

133.9 

351.0 

2.6214 

1.619 

1.61 

1.80 

Na  Br 

102.9 

366.1 

3.5568 

1.886 

1.81 

2.03 

NaCl 

58.5 

410.0 

7.0085 

2.647 

2.41 

[40] 

2.70 

2.38 

Nal 

149.9 

297.0 

1.9813 

1.408 

1.51 

1.69 

RbBr 

165.4 

385.0 

2.3277 

1.526 

1.39 

1.56 

Rb  Cl 

120.9 

444.0 

3.6725 

1.916 

1.66 

1.86 

RbF 

104.5 

490.0 

4.6890 

2.165 

2.13 

2.39 

Rbl 

212.4 

331.0 

1.5584 

1.248 

1.17 

1.31 

Iron 

55.9 

415.3 

[52] 

7.429 

2.726 

3.85 

[44] 

4.31 

3.16 

PMMA 

97.1 

bo 

2.2842 

1.51 

1.373 

[1] 

1.54 

1.36 

Uc  data  from  Reference  53,  except  as  noted. 
Ct  data  from  Reference  52,  except  as  noted. 
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TABLE  B-3.  Comparison  of  Elastic  Wave  Velocities  and  Heat  of  Detonation 
for  Selected  Energetic  Materials. 


Substance 

AHdet 

MJ/KG 

AHdet 

(km/sec)^ 

V  AHdet 

km/sec 

Cb 

km/sec 

Cs 

km/sec 

Comp-B3 

4.69(E) 

4.69 

2.166 

2.42 

1.612 

Comp-B 

5.86(C) 

5.86 

2.421 

2.42 

1.612 

PETN 

5.73(E) 

5.73 

2.394 

2.32 

1.582 

PBX-9407 

6.11(C) 

6.11 

2.472 

2.32 

1.070 

PBX-9404 

5.36(E) 

5.36 

2.315 

2.26 

1.491 

TNT(P) 

4.27(E) 

4.27 

2.066 

2.08 

1.310 

TETRYL 

4.56(E) 

4.56 

2.135 

1.76 

1.165 

TATB 

4.52(C) 

4.52 

2.126 

1.44 

1.000 

Octol 

5.98(C) 

5.98 

2.445 

2.49 

1.601 

RDX 

5.94(E) 

5.94 

2.437 

2.65 

PBX-9501 

6.03(C) 

6.03 

2.455 

2.50 

~ 

PBX-9011 

5.69(C) 

5.69 

2.385 

2.41 

- 

PBX-9010 

5.69(C) 

5.69 

2.385 

2.13 

~ 

LX-I7 

4.27(C) 

4.27 

2.066 

2.04 

(P)  -  Pressed 
(C)  -  Calculated 
(E)  -  Experimental 
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Figiue  B-1.  AHdet  or 


AhpET 

niAv 


Versus  Cb- 
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CALCIUM  ATOM  EFFECTS  ON  BONE 
VIBRATION  RESONANT  FREQUENCIES 

In  certain  bio-mechanics  investigations  with  vibrating  bone  segments  [60],  it  was  found 
that  the  frequency  relation  given  by  Fitzgerald  [2,  Chapter  II]  was  the  one  that  fit  the  data.  This 
frequency  relation  is: 


Vm  = 


(0-1) 


where 


h  =  Planck’s  Constant 

m  =  Mass  of  the  atom  responsible  for  the  resonance  frequency 
s  =  Repetitive  segment  length  of  the  structure 
n  =  integer,  n  =  1, 2, 3,  etc. 

The  resonant  frequencies  were  sensitive  to  the  mass  of  the  calcium  atoms.  In  fact,  Fitzger¬ 
ald  [61]  pointed  out  that  the  mass  for  an  isotope  of  calcium  provided  a  better  comparison  with 
the  data.  This  observation  w  is  substantiated  by  the  bio-mechanics/medical  researchers  [62]. 

If  the  influence  of  particular  atoms  appears  in  such  a  sensitive  manner  in  vibrating  bone,  it 
is  to  be  expected  that  the  sensitive  influence  of  certain  atomic  elements  will  be  found  in  shock- 
loaded  materials  such  as  polymers,  compounds,  and  mixtures  of  compounds. 

If  Equation  C-1  is  multiplied  by  h,  it  is  equivalent  to  Equation  12.59  in  Reference  41,  for 
the  energy  of  a  particle  in  a  one  dimensional  box  (or  line  segment). 
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